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Super-heavy Quarkonia as a probe of the QCD Vacuum 

J. Fingberg, University of Wuppertal, 42097 Wuppertal, Germany 

We investigate the effect of a non-perturbative gluon condensate on the spectrum of bound states of extremely 
heavy quarks using Non-Relativistic-QCD (NRQCD). Simulations are done with an improved gluon action on 
small lattices. We compare our numerical results with perturbative calculations to test the consistency. 



The phenomenon of gluon condensation is un- 
der vital discussion since more than fifteen years 
now jjj. The gluon condensate is an important 
quantity not only in the sum rule approach but 
also in high-energy hadron scattering and phe- 
nomenological models of confinement. Despite 
many indications of a non-trivial vacuum the 
quantitative picture is still vague ||]. Quarkonia 
are the simplest hadronic systems that allow to 
test the vacuum structure of QCD. Heavy quarks 
can be treated as non-relativistic objects in a glu- 
onic vacuum making them tangible to numerical 
simulations and perturbative calculations. 



In this study we compare the quark mass depen- 
dence of quarkonium level splittings from sim- 
ulations using NRQCD 0] to analytic results. 
The effect of the gluon condensate can be calcu- 
lated perturbatively assuming that super-heavy 
quarkonia are almost Coulombic: E n i = 2mq + 
E^ai-mer _|_ _ A non- vanishing gluon conden- 
sate § (E 2 ) = -1/4(0|G^G^|0) + gives rise 
to a level shift SE n i by quadratic Stark split- 
ting. The simplest ansatz is a constant chromo- 
electric field (|||. ^ n this case the shift de- 
pends on the principal quantum number and the 
quark mass in the form 5E n i oc n 6 /mg. A fi- 
nite gluon correlation length [||J^] (Ac? — 0.22 fm 
weakens this dependence: 5E n i oc n 4 /mg. 
There are two length scales in the system, Aq, 
and the quark correlation length, Aq, thus two 
limiting cases can be considered jl2|: the non- 
potential case Aq >> Aq and the potential case 



Ag << Aq . The first case reduces to Leutwyler's 
result, AE 2 i = E 2 q-E w = ciitlq / n 2 + c 2 n 6 / 
as uiq — > oo. For the second case one finds 
AE 2 i — c^mq/n 2 + C47t, 4 /toq. The hyperfmc 
splitting of the ground state can be parametrized 
in the same form. 

Our simulations are done with a Symanzik and 
tadpole improved gauge action with a plaqucttc 
and a rectangle term 
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enabling us to use a small lattice of size 16 4 . 
The quark propagators are computed with the 
NRQCD evolution equation: 



G t+ i = (1 - 



2n 



2n 



n (l-a8H) G t . 



The Hamiltonian includes all spin-independent 
interactions to 0(v 2 ) and spin-dependent interac- 
tions to 0{v ). Difference operators are improved 
to avoid large discretization errors. We generate 
« 300 decorrelated configurations at f3 = 9.17 
and determine the meson spectrum for S and P 
states. To set the scale we measure the heavy 
quark potential for on axis and \J~2 separations 
from Wilson loops and determine the lattice spac- 
ing from the force [@, F(R) = dV(R)/dR. A 
comparison of dimensionless quantities 7q F (ro) 
on the lattice and using both Eichten's phe- 
nomenological potential, 



V(r) = -n/r + err , y/a = 427 MeV, 



0.52 



and Liischer's bosonic string model, — tt/12, 
shows good agreement of the sting model param- 
eterization with the lattice data. A fit shown in 
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(Fig. 1) 

r_0 = 0.33 fm 

1/a = 4.6 +/-0.1 GeV 
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fig.l gives the results tq = 0.34 fm and a -1 = 
4.6(1) GeV. Meson propagators and effective 
masses are measured for 4 values of the bare 
quark mass, rriQa = 1.5,2,3,4, using Coulomb 
wave functions. The kinetic mass is extracted 
from the dispersion relation, m(p) = m s i m + 
\p \ 2 /2rrikin — \p | 4 /8m2, by a fit to 1 S finite 
momentum meson masses. In fig. 2 we see that 
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(Fig. 2) 
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for low rriQ the results nicely extrapolate to the 
cc and bb points of The splitting AE21 = 

E2( l S )-Ei( 1 Sq) plotted in fig.3 is determined 
from correlated multi-exponential fits. The ex- 
perimental values also indicated in fig.3 are 670 
McV for cc and 558 Mev for bb. The important 
observation is that the new data points continue 
to decrease with increasing quark mass. The ef- 
fect becomes smaller with increasing quark mass 
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and vanishes in the noise for the two largest 
masses. The hyperfine splitting can be deter- 
mined very accurately from a fit to the ratio of 
3 S\ and 1 Sq meson propagators with a single ex- 
ponential. The experimental value for AE^ S is 
118 MeV. We observe a steep decrease with quark 
mass which levels off at higher values of rriQ. 
From analytic calculations we expect a behaviour 
of both splittings, AE21 and AE^fs, proportional 
to the quark mass for a pure Coulomb potential. 
A gluon condensate would modify this result to 
AE = ciiriQ + C2(a s G 2 }/mQ, with an exponent 
a = 2 ... 3 depending on the relative magnitude 
of Ag and Xq. In fig. 5 we see qualitative agree- 
ment between our data and analytic predictions if 
the value of the gluon condensate is large enough. 
The observed decrease of the splittings does not 
agree with the expectation from a pure Coulomb 
potential. The apparent contradiction is likely to 
prevail for any reasonable potential model calcu- 
lation p3|. Finally the ratio of the 1 Sq and the 
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hyperfine splitting shows an increase at least up 
to m,Q = 15 GeV while a pure Coulomb potential 
would give a constant. The experimental value is 
r ce = AE 21 /AE hfs = 5.7. 

SUMMARY 

We present results of a NRQCD analysis based 
on simulations with a Symanzik and tadpole im- 
proved gauge action on a small lattice. We find 
that E2( 1 So)—Ei( 1 Sq) and the hyperfine splitting 
of the ground state decrease with increasing tuq 
in accord to existing data for cc and bb systems at 
small vtiq. The decrease is in qualitative agree- 
ment with model calculations of gluon condensa- 
tion if the gluon condensate is large enough. How- 
ever, present systems are still far from Coulombic. 
A direct determination of (a s G 2 ) is difficult and 
we need even higher quark masses itiq > 20 GeV 
and also higher excitations n — 3,4 (in prepara- 
tion). A promising signal which requires extreme 
precision is the S-P splitting and P-state hyper- 



fine splitting as noted in |14|]. On the analytical 
side an improved perturbative calculation of the 
corresponding matrix elements is desirable. 
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